Kondo lattice model with a direct exchange interaction between localized moments 
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We study the Kondo lattice model with a direct antiferromagnetic exchange interaction between 
localized moments. Ferromagnetically long-range ordered state coexisting with the Kondo screening 
shows a continuous quantum phase transition to the Kondo singlet state. We obtain the value of 
the critical point where the magnetizations of the localized moments and the conduction electrons 
vanish. The magnetization curves yield a universal critical exponent independent of the filling fac- 
tors and the strength of the interaction between localized moments. It is shown that the direct 
exchange interaction between localized moments introduces another phase transition from an anti- 
ferromagnetic ordering to a ferromagnetic ordering for small Kondo exchange interaction. We also 
explain the local minimum of the Kondo temperature in recent experiments. 
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I. INTRODUCTION 

Interaction between itinerant electrons and localized 
moments in the Kondo regime is described by the s- 
f mixing and the polarization terms in the Kondo lat- 
tice model. The polarization term induces a RKKY type 
magnetic correlation, while the s-f mixing causes demag- 
netization. The magnetic correlation is dependent on the 
filling factor; ferromagnetic for low filling and antiferro- 
magnetic around half filling. 1 As the interaction between 
the conduction electrons and the localized moments in- 
creases, the phase evolves to the Kondo singlet phase. 
In particular, for half-filled case, the phase transition has 
been shown continuous and the value of critical point has 
been obtainedi 2 i 3 i 4 i 5 

In reality, however, the direct exchange interactions be- 
tween localized moments should be included for a more 
comprehensive understanding. The Kondo exchange in- 
teraction, Jk , between the conduction electrons and the 
localized moments invokes the effective exchange inter- 
action between the localized moments, the RKKY in- 
teraction Jrkky ~ \Jk\ 2 ' 6]7]8 In the present study, in 
addition to the RKKY interaction, we introduce a direct 
antiferromagnetic exchange interaction between localized 
moments. In the ferromagnetic Kondo lattice model for 
manganites, its effects have been widely investigated^ In 
this study, we aim to address the effect of antiferromag- 
netic exchange interaction, J a f , on the phase diagram and 
on the critical behavior of Kondo lattice model (KLM) 
for low filling cases by calculating the density of states in 
a mean field level. 

First of all, we investigate the magnetizations of both 
the lattice magnetic moments // and the conduction elec- 
trons fj, for various filling factors and various values of J a f . 
The magnetic moments //and fj. show qualitatively differ- 
ent behaviors; [t! falls monotonously from its maximum 
value as Jk increases, whereas fx starts from zero, passes 
through a local maximum, and falls to zero at the same 
point as \x' does. Near the critical point, the magnetiza- 
tion curves show characteristic behaviors of continuous 
phase transition. The critical behavior of the magneti- 



zations shows a universal exponent independent of the 
filling factor and the strength of exchange interaction J a f 
such that fj/(fj.) ~ {Jk — Jc) 1 ^ 2 with the critical point 
J c . There emerge two mean fields. The first one, x, rep- 
resents the mean field of the s-f mixing and the other 
field, B, the magnetic correlation between localized mo- 
ments. As Jk increases, the x field increases whereas 
the B field decreases, which indicates a phase transition 
from the magnetically ordered state to the Kondo state. 
Calculated values of the mean fields x and B explain con- 
sistently the the physical situation such as the existence 
of phase transition between the ferromagnetic ordered 
state and the Kondo singlet state. 

A notable feature of the present calculation is the ex- 
istence of a new antiferromagnetic state at small Kondo 
exchange interaction especially for low filling factors. 
The s-f mixing field, x, is generally believed to be a 
monotonously increasing function of Jk^& However, for 
very small Jk, we find a local minimum of x field. It 
is accompanied by suppression of the B field, which is 
the precursor of a transition to another magnetic order, 
an antiferromagnetic phase. A recent experiments on 
CeRu2Ge2 by Siillow et ali& and on CeAgSb2 by Sidorov 
et ali^ can be explained in connection with the anoma- 
lous behavior of the x field and the B field. 



II. COEXISTENCE OF MAGNETIC ORDERING 
AND KONDO SCREENING 

The magnetic correlation and the s-f mixing process 
are the two competing orders in the KLM. The Doniach 
diagram 12 shows that the former has dominant effect for 
small Jk while the latter for large Jk- The Kondo ex- 
change model includes the s-f mixing term and the po- 
larization term in the Kondo interaction to describe the 
competing orders. The polarization term give rise to the 
RKKY interaction^ 

Integrating out the fast mode in the action including 
the Kondo interaction may produce the effective RKKY 
interaction term in the Hamiltoniar_ 13 i with the interac- 
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tion strength dependent on the Kondo interaction and 
the cutoff frequency. In this study, however, we con- 
sider the case that the RKKY interaction is implicitly 
included in the polarization term of the Kondo Hamil- 
tonian. Instead we introduce the direct antiferromag- 
netic exchange interaction between the localized mo- 
ments whose strength we treat as independent parameter. 
For ferromagnetic Kondo lattice model for manganites, 
where the interaction between the conduction electrons 
and the localized electrons is ferromagnetic, the effect 
of the direct exchange interaction between localized mo- 
ments has already been studied^ Here, various antifcrro- 
magnetic phases appear, even though the antiferromag- 
netic exchange interaction is much weaker than that of 
the Hund coupling 

The Hamiltonian H = Hq + Vk + Vh with the Kondo 
interaction Vk and the interaction between localized mo- 
ments Vh is written by 



moments respectively&iSii& Hereafter, we consider the 
uniform case such that x\ = x and y\ = y. 

Here, H± describes the Kondo scattering and H2 the 
polarization of the conduction electrons by the localized 
moments. H2 gives the magnetic ordering through the 
RKKY interaction between localized moments. The com- 
petition between H\ and H2 is previously studied^. H3 
comes from the transverse part of the direct exchange 
interaction Vh in Eq. JSJ and can be represented in k- 
space^ 



H 3 = Bj2^ftfr«, 



(8) 



k.a 



Ho = 



(i) 



Vk = ~^E ( Si ' Si - | n /4 n c,ij ( 2 ) 

Vh = ~~Jai E ( Si ' Si + a ~ 7 n /.i n /.i+a ) . (3) 



where o represents spin, a a unit vector connecting the 
nearest neighbors, and Eq the energy of the localized 
level. Jk(< 0) is the antiferromagnetic Kondo exchange 
interaction and J a f (< 0) the antiferromagnetic exchange 
interaction between localized moments. ct CT is the cre- 
ation operator for conduction electrons with wave vector 
k, ]\ a for localized electrons at cite i, n c i — J2a c to- c icn 
and n fA = £ a _ 

Using the Stratonovitch-Hubbard transformation and 
the identities, Si = §C;<TCj and Si = |/j c/i, with Pauli 
matrices a, the interaction terms are represented by Gft-(ui) = N , 

V K = H 1 +H 2 - (J K /2) x? and V H = H 3 + H 4 - ^ 

(Jaf/4)Ei,aJ/. ? . Where 

H > = 4rE-'-C'- • '•• ,.A. -> (4) 

i,cr 

H * = ^E 4 ,h. ■ 4,^ • 4 „ G cT M = ^ 



where B = —zyJ a f/2D and ek = — (D/z) ^ a cosk • a 
with coordination number z. This term gives the local- 
ized moments a small band width 2BD. if 4 is the term 
polarizing the nearest neighbor localized moments. 

In order to describe the critical behavior near the crit- 
ical point, we consider the coexistence of the magnetic 
ordering and the Kondo screening. Calculating magne- 
tizations of the localized moments and the conduction 
electrons for various filling factors and various J a f, we 
obtain the phase boundaries and the critical exponent. 
In this case we have the relations such that 

( c i| c U> = \( n + M), <4 T c iT> = \{n - M), (9) 

(/i/iiH^l-A (/ I t T /lT) = 5(l+M'). (10) 

where /i (//) is the magnetization of s (f) electrons in the 
spin-down (spin-up) direction and n is the filing factor. 
Since the Kondo exchange interaction Jk is antiferro- 
magnetic, we expect that [i and // have the same sign. 

With these equations the Green's function for the f- 
electrons and the c-electrons corresponding to the Hamil- 
tonian H can be obtained as follows, 



oj ± is — Eq — Be k - -JK(n + (i)— 



72 2 



-zJ af {l - fjf) - . 
4 4[w ±is - e k - - M )J 



(11) 



(5) 
(6) 



u> ± is - e k - -Jjc(1 - M')" 



J K x 2 



4[w ± is - £Jo - Be fc - j Jif(n + m) - |zJ af (1 - //)] 



(12) 



J, 



af 



^ ] [(i i.crfuo) 4fa,— (7 /i+a.-er 



For the spin down case, the quantities such as G/|(w) 
can be simply obtained by substituting /i — fj, and 



/j t CT fi , a ( /}_(_ a t _ o- /i + a , - ct ) 
(/i,cr/ i . cr ) (/i+a.-CT/i+a.- ")]' 



(7) 



The above Green's functions can be arranged to be 
written as follows, 



where xi and yi are the mean fields representing the s-f 
mixing and the magnetic correlation between localized 



f// T (w) ± is - e k hfi(u)±is- e k 



3 



(13) 



g cT (w) ± is - efe /i cT ± is — e k 



(14) 



Here we use the flat band for the conduction elec- 
trons such that Im^ fc l/(cj ± is — e k ) = 1/25(0) for 
—D < oj < D (otherwise). Hence, the density of 
states, pfi(uj) = ImG/|(oj), can be obtained by pf-\(oj) = 
5/ T (w)/25 (Qf}(bj)/2D) for w/if < u < W/2T (^/3T < 
uj < lo ) . The density of states for the conduction elec- 
trons, p c ^(uj) = ImG c |(o;), can be obtained in a similar 
fashion. Since we can readily show Pf^(u>) = Qf^(u>) 
and P c t( w ) = Qc]{u) in the above, the density of states, 
P/t( w ) = /5/t(^t) and /MM = Pct(^t)> are S iven b y 



p/t(^t) 
Pct(^t) 



455 
1 

45 



o 



02 + BJ 2 x 2 



o. 



O 2 + BJ 2 x 2 



(15) 
(16) 



where 



fi T (w) = (l-B)uj--J((n + p)-(l-p')B) 

- 5 0T - \zJ &i (1 - //). (17) 

The band edges can be obtained by solving the equa- 
tions — D < gf( c ^(uj) < D and — D < hf( c ^(uj) < D. 
Using the expressions obtained for g^(w) and h^(ui), 



= 7Td 



1 

25 



(l+B)LU--J((n + p) + (l-p')B) 



5, 



ot 



?Jaf (1 - M') 



O 2 + J 2 2; 2 5 



(18) 



(1 + B)w--J((n + ^ + (1-^)5) 



5, 



izJ af (1 - p!) - ^Jn 2 + J 2 x 2 B 



,(19) 



we get for u>^ 



1 



wit 



W2| 



W3| 



W 4 f 



[u a -(l + B)D 



2 

V {oJb 
1, 

l r 



f- (1 - B)D} 2 - ABD 2 + J 2 x 2 ] (20) 
(1 + B)D 

- (1 - 5)5} 2 - ABD 2 + J 2 x 2 } (21) 
(1 + B)D 

b (1 - 5)5} 2 - 455 2 + J 2 x 2 ] (22) 
(1 + 5)5 

- (1 - 5)5} 2 - 455 2 + J 2 x 2 ], (23) 



where w Q = \.J{(n + p) + (1 -//)) + #c , .. 
and uj b = \ j{{n + p) - (1 - //)) + 5 0T + ±zj af (1 - p'). 
Here, we used the relation, gf a (uj) = g ca (uj) = g a {u), 
hfo-fo) = h ca (uj) = ha-(u), and Uf^ = w c i\ = <^i\- 



re w a = jJ((n + p) + (1 -/i 

Wfc - \ j{{n + n) - (1 -p')) 
3. we used the relation. Of, 



-0| + jZJa£(l _ /"') 



III. RESULTS AND DISCUSSIONS 

The ferromagnetic magnetizations p' and p in presence 
of the Kondo screening can be obtained as a function of 
the Kondo exchange interaction Jk by solving the self 
consistent equations, Eqs. © and 110|) . using the density 
of state, Eqs. |(TSJ and (JTSJl. Eqs. © and (JTDJl for spin 
up case are represented as 



1 + // 



n — p 



1 - B 



1 -5 



/5 c j(0)<iO. 



(24) 
(25) 



From the above equations Oft an d ^it can be evaluated 
by inserting the density of states from Eqs. (|15fl and 
H16| l. Oj(5f) and 0-|-(o;it) can a ^ so be represented from 
Eq. (fTTjl . Solving the equations, Oft = 0|(5f) and 
fii-f = n-f^i-j-), -Eqt an d 5ft can be gi yen by 



40 2 T 



{I- B) 2 {ABD 2 - J 2 x 2 ^ 



2(1 + 5)011 + (1 - B) A /40 2 - 16B 2 D 2 + ABJ 2 x 2 



T 

-i J((n + M) - (1 - MO) - J^af (1 - M') - (1 - 5)5 

(26) 



5ft — 



1 - 5 



fO 



ft 



5 



01 



+ 1 J((n + m) - (1 - A*')^) + j^af (1 - M') 



(27) 



The values of magnetizations, /i' and ^t, are obtained 
numerically. First, we fix the values of x and B for 
given values of n, J a f , and Jr- and, then, solve the equa- 
tions, 5ft = Ef\ and 5 



5ox, to obtain the values 
of magnetizations, /x' and p. From these values we cal- 
culate the total energy of the conduction electrons and 
localized moments. The total energy can be written as 

5 = Yjc Sui a ^/"-M + Pca(^)]duj + E M f- E M f can 
be obtained from the Hamiltonians in Eqs. I0J-0 with 
the mean values from Eqs. © and l|l()(l and is given by 
Emf = -(1/2) J K x 2 - BD 2 /zJ af - (1/4) J K {n + pp 1 ) - 
(l/8)zJ a f(l — p' 2 ). For appropriate ranges of x and 5 
values, we repeat the above process and determine the 
values of p! and p at the minimum point of the total 
energy of the system. 

In Fig. ^ we show the magnetization of the localized 
moments, p! , and of the conduction electrons, p, in fer- 
romagnetic state when J a f/5 = —0.001. The magnetiza- 
tion curves clearly demonstrate existence of continuous 
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FIG. 1: Magnetization of the localized moments, \j! and of the 
conduction electrons, fi in ferromagnetic state for various fill- 
ing factors, n, as functions of the Kondo exchange interaction 
strength, Jk- Here, we set the strength of direct exchange 
interaction J a f/D = —0.001. 
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FIG. 3: Log- log plot of the magnetization of localized mo- 
ments fi' and conduction electrons fi near the critical point 

Jo- 



quantum phase transition. We observed that the criti- 
cal values of the Kondo exchange interaction are higher 
for smaller values of filling factors, which is consistent 
with previous studies^ However, there have been con- 
flicting reports on the behavior of fi. A recent Monte 
Carlo study^ showed existence of a local maximum in the 
/i curve, but its behavior for small values of Jk was left 
uncertain. On the contrary, a more recent bond-operator 
mean-field study^ produced a monotonously decreasing 
behavior. Present study shows that the magnetization of 
conduction electrons \i starts to increase from zero and 
pass a local maximum before it falls to zero at the same 
points as those of fjf. It is consistent with the Monte 
Carlo calculation and also clarifies the situation for small 
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FIG. 2: Magnetization of the localized moments, // for vari- 
ous values of direct exchange interaction J a f between localized 
moments with a filling factor, n = 0.2. 



values of Jk- 

Fig. [21 shows the behavior of \i' near the critical point 
for n=0.2 with various values of J a f. J a f competes with 
the ferromagnetic correlation and, thus, makes the criti- 
cal points slightly smaller as shown in Fig. [21 We observe 
that the critical points are dependent on the values of J a f 
as well as the filling factors. Close study on the critical 
behaviors of n' and fj, clearly shows that the critical ex- 
ponent has a universal value independent of the filling 
factor and J a f. The form is given by (Jk — Jc) 1 ^ 2 as 
shown in Fig. [21 

The mean field variable x representing the s-f mix- 
ing is closely related to the Kondo temperature and the 
variable B representing the correlation between localized 
spins is to the correlation temperature^ The behaviors 
of x and B as functions of Jk are plotted in Fig. Ufa). 
Fig. 01a) shows that both x and B have higher values for 
higher filling factors, which is consistent with the study 
by Ruppenthal et alii The field B is shown to behave in 
the opposite direction of x. This signifies that B repre- 
sents the magnetically ordered state competing with the 
Kondo singlet state. 

The KLM with a direct exchange interaction between 
localized moments exhibits a qualitatively different be- 
haviors from those without direct interaction between 
localized moments. The value of x shows a local min- 
imum as Jk becomes quite small. This phenomenon is 
accompanied by the suppression of the B field. This be- 
havior indicates that the magnetic correlation becomes 
weak and another transition is being anticipated. The 
local minimum of x is shown in Fig. 01 (a) for various fill- 
ing factors when J a f/D = —0.001 and for various values 
of J a f when rt=0.35 in (b), where the re-growth of x is 
shown stronger for larger |J a f |. 

We have calculated the antiferromagnetic state energy 
and compared it with that of the ferromagnetic state. In 
Fig. [SJ the phase boundary for the ferromagnetic state 
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and the antiferromagnetic state is shown in the lower part 
of the figure. For low filling, the antiferromagnetic state 
become more favorable. As the value of Jk increases, the 
transition to the Kondo singlet state follows as shown 
in Fig. [5] Similar phase boundary has been obtained 
previously by comparing the energy of the Kondo singlet 
state with that of the ferromagnetic state without Kondo 
scattering. 1 Here, it should be noted that the present 
calculation is carried out when the ferromagnetism and 
the Kondo screening coexist. 

This behavior explains the recent experiments on 
CeRu2Ge2 by Siillow et aliifi and on CeAgSb2 by Sidorov 
et a/^i In these experiments, the interaction between the 
localized moments, Ce, is ferromagnetic while the RKKY 
interaction is antiferromagnetic. Two consecutive phase 
transitions have been observed: a transition from the fer- 
romagnetic phase to the antiferromagnetic phase followed 
by a transition to the Kondo singlet state as Jk increases 
further. If we exchange the role of the antiferromagnetic 
interaction with that of the ferromagnetic interaction, we 
can explain the experiment in the present scenario, since 
here we have considered the antiferromagnetic interac- 
tion between f-electrons and the ferromagnetic RKKY 
interaction. 

In the experiment on CeRu2Ge2^ we observe suppres- 
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FIG. 4: Plots for the mean field x and B (a) for various n 
values when J a f/D — —0.001 and (b) for various values of 
Jaf/D when n = 0.35. 
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FIG. 5: Phase diagram for various J a f. The phase boundary 
for the antiferromagnetic state and the ferromagnetic state is 
shown in the lower part of the figure. The phase boundary for 
the ferromagnetic state and the Kondo singlet state is shown 
in the upper part for Jat/D — and J^t/D = —0.001. 



sion of the magnetic transition temperature near the crit- 
ical point where the transition from one magnetic state to 
another magnetic state occurs. This corresponds to the 
suppression of the B field in the present study. We can 
see in the experiments the local minimum of the Kondo 
temperature, Tk- In the experiment on CeRu2Ge2^ we 
can only observe that the Kondo temperature resists go- 
ing to zero as Jk decreases. However, in the experi- 
ments on CeAgSb2»ii and on CeRhlngAi, we clearly see 
the local minimum of the Kondo temperature. This local 
minimum has been conjectured to be caused by a cer- 
tain mechanism competing with the Kondo effect. In the 
present theory, we show that the direct interaction be- 
tween the localized moments competing with the Kondo 
scattering gives rise to the local minimum of the x field 
and, thus, the Kondo temperature. 



IV. CONCLUSION 

We have studied the KLM with an antiferromagnetic 
exchange interaction between localized moments for a 
general situation in which ferromagnetic correlation co- 
exists with the Kondo screening. Calculating the mag- 
netization of the localized moments and the conduction 
electrons as a function of the Kondo interaction for var- 
ious filling factors and values of J a f , we obtained a con- 
tinuous quantum phase transition from the ferromagnetic 
ordered state to the Kondo singlet state and the critical 
points where the magnetizations vanishes. Around the 
critical points, we obtained the critical exponent inde- 
pendent of the filling factor and the values of J a f such 
that m' (m) ~ {Jk ~ Jc) 1/2 - 

Furthermore, we obtained the behavior of the mean 
fields, x and B, corresponding to the Kondo scatter- 
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ing and the magnetic correlation respectively. As the 
Kondo interaction increases the value of B decreases 
while x increases, which explains existence of the Kondo 
singlet state at large value of Jk- On the contrary, 
as Jk decreases, the value of B increases while x de- 
creases. Hence, the enhanced magnetic correlation with 
suppressed Kondo screening leads the magnetically or- 
dered state. 

When Jk is further reduced, the another phase tran- 
sition from the antifcrromagnetically ordered state to 
the ferromagnetically appears due to the J a f interaction. 
Comparing the energy of the ferromagnetic state and the 
antiferromagnetic state, we obtain the phase boundary, 



which shows clearly that an antiferromagnetic state ap- 
pears for low filling. The x field exhibits a re-growth ac- 
companied by the suppression of B. This phenomena is 
due to the J a f interaction and is a precursor for the tran- 
sition to the antifcrromagnetically ordered state. Recent 
experimentsi2iii*i£ can be explained through the present 
theory. 
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